Does the neutrino magnetic moment have an impact on solar neutrino physics? 
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Solar neutrino observations coupled with the recent KamLAND data suggest that spin-flavor 
precession scenario does not play a major role in neutrino propagation in the solar matter. We 
provide approximate analytical formulas and numerical results to estimate the contribution of the 
spin-flavor precession, if any, to the electron neutrino survival probability when the magnetic moment 
and magnetic fleld combination is small. 
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I. INTRODUCTION 

In the earlier days of the solar neutrino research ac- 
tivities one of the more speculative solutions proposed to 
resolve the puzzle of missing neutrinos invoked the inter- 
action of the neutrino magnetic moment with the solar 
magnetic fields. Although initial attempts H ignored 
matter effects, eventually the combined effect of matter 
and magnetic fields was brought out 0, Q ■ Simultaneous 
presence of a large neutrino magnetic moment, magnetic 
field combination and neutrino flavor mixing can give rise 
to two additional resonances besides the MSW resonance 
^1^. Initial numerical calculations using the resonant 
spin-flavor precession scheme were carried out before the 
gallium experiments were completed. These calculations 
hinted a solution of the solar neutrino problem with the 
neutrino parameters in the LOW region provided that 
there is transition magnetic moment as large as 10~^^^b 
and a magnetic field of the order of 10^ G. With this solu- 
tion count rate at the gallium detectors would be signif- 
icantly reduced. Even though gallium experiments ruled 
out this particular solution variants of the spin-flavor pre- 
cession solution to the solar neutrino problem continue to 
be investigated by many researchers. (A representative 
set of the recent work is given in Refs. |3jIEHiIi3)- 

In the meantime experimental data has increasingly 
disfavored the spin-flavor precession solution to the so- 
lar neutrino problem. Earlier reports of an anticorrela- 
tion between the solar magnetic activity and solar neu- 
trino capture rate at the Homestake detector (lli] was a 
prime motivation for considering the magnetic field ef- 
fects. An analysis of the Su per -Kamiokande data rules 
out such an anticorrelation |l2j |. (For a counter argu- 
ment, however see Ref. Also, if the neutrinos are 



t Permanent Address 
'Electronic address: 
■f Electronic address: 



baha@nucth.physlcs.wisc.edu| 
volpe@ipno.in2p3.tr| 



of Majorana type, such a scenario would produce solar 
antineutrinos 0, 0| . Both the Super-Kamiokande 
and KamLAND collaborations [T^ failed to find any ev- 
idence for solar antineutrinos. Finally a global fit of all 
solar neutrino experiments without using the spin-flavor 
precession solution (see e.g. Ref. '^^s confirmed by 
the KamLAND experiment. 

On the other hand our knowledge of both the solar 
magnetic fields and neutrino magnetic moments has been 
improving. (At the very least we now know that that neu- 
trinos, since they are definitely massive, have magnetic 
moments the magnitude of which depends on the physics 
beyond the Standard Model). Hence it may be worth- 
while to revisit the spin-flavor precession mechanism and 
knowing it is not the dominant mechanism, explore what 
implications it may still have. 



II. STATUS ON SOLAR MAGNETIC FIELDS 
AND fi„ 

Neither the magnetic pressure in the core nor the im- 
pact of the structure of the magnetic field on the stellar 
dynamics are usually taken into account in the Standard 
Solar Model 0, l20| . Within the uncertainties of the nu- 
clear physics input [2j Standard Solar Model agrees well 
not only with the neutrino observations but also with the 
helioseismological observations of the sound speed profile 
Hit 1^ • Direct measurements of the magnetic fields 
in the radiative zone with acoustic modes are not possi- 
ble, however even with a one-dimensional solar model a 
sizable magnetic field would contribute additional pres- 
sure. An analysis found that a magnetic field greater 
than 10^ G localized at about 0.2 Kq would cause the 
sound speed profile to deviate from the observed values 
[2^. However this particular calculation does not place 
restrictions on the magnetic field exactly at the center 
of the Sun or even at O.I Kq. Similarly magnetic field 
strengths greater than ~ 7 x 10^ G are not allowed since 
they could cancel the observed oblateness of the Sun, 



2 



making it spherical and even prolate j25l| . 

Helioseismology provides more detailed information 
about the magnetic fields in the convective zone. The 
rotation profile of the Sun is presently known down to 
about 0.2 R0: the rotation of the solar radiative zone 
is like that of a solid body rotating at a constant rate 
[26l|. Such a rotation profile suggests the existence of a 
magnetic field in the radiative zone. The observations of 
the splittings of the solar oscillation frequencies can be 
used to infer the magnetic field. The odd terms in the 
azimuthal order are determined only by the rotation rate 
in the solar interior. The even terms may receive addi- 
tional contributions from the magnetic fields. An analy- 
sis of the helioseismic data indicates that rotation alone 
is not sufficient to explain the observed even splitting co- 
efficients. Other helioseismic observations are consistent 
with a magnetic field of ~ 20 G at a depth of 30000 km 
below the solar surface Similar arguments limit the 
toroidal magnetic field to < 300 kG at the bottom of the 
convective zone. On the contrary, at present, there is no 
direct helioseismic evidence for the presence or absence 
of sizable magnetic fields in the radiative zone. 

As far as the neutrino magnetic moment is concerned, 
best present direct upper limits come from reactor ex- 

timents, i.e. /i^ < 1.0 - 1.3 x IQ-^^fXB at 90% C.L. 
[3^ which improve previous bounds |^ l3^ 133 . l33| , 
as well as from the recent Super-Kamiokande solar data 
which has given the limit of < 1.5 x 10^^"/i_B at 90% C.L. 
|35j| . Combining recent solar neutrino experiments with 
the KamLAND data yields a limit of < 1.1 x 10^^"^b 
at 90% C.L.. From astrophysical considerations, an in- 
direct upper limit in the range of 10~^^ — 10~^^/iB have 
been obtained j37j . the exact limits being model depen- 
dent. New experiments are now under study which would 
lower the direct limits down to the level of a few x 10""'^^ 
using a static tritium source [sl, |3^ ^3 while the use of 
low energy beta-beams might lower it by about one order 
of magnitude |40|, |4l| . 



III. SPIN-FLAVOR PRECESSION FORMALISM 
AND RESULTS 

The evolution of the chiral components of two fiavors 
of neutrinos is described by Q 



(3) 



down for the left-handed component 

hW=( ^sin2^? \ 

For the Majorana neutrinos the right-handed part of the 
Hamiltonian, is given by replacing Ve and in Eq. 

© hy —Ve and —Vfj, respectively. In these equations the 
matter potentials are 



Gp 
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and 



(2iVe-iV„), 



(4) 
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where Gf is the Fermi constant of the weak interactions, 
Nf, is the electron density, and iV„ is the neutron density. 
In the above equations for the Dirac neutrinos a general 
magnetic moment matrix is possible: 



M 



f-^ee l-^efi 



(6) 



For the Majorana neutrinos the diagonal components of 
Eq. (|SJ) vanish and the off-diagonal components are re- 
lated by -fiefi = Mmg = 

In this scenario there are several resonances. In addi- 



tion to the standard MSW resonance 
takes place where the condition 

r 2 

V2GFNe = — COS 261 



) that 
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is satisfied in the Sun. For the left-handed electron neu- 
trinos that are produced the core of the Sun a second. 



spin-flavor precession, resonance (i', 



) is pos- 



sible. For the Dirac neutrinos it takes place where the 
condition 



Gf 
V2 
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is satisfied whereas for Majorana neutrinos it is where 
the condition 



< T\ 



dt 



(L) 



For the Dirac neutrinos one has 
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and H^-^-^ is given by setting and equal to zero in 
Eq. For the Majorana neutrinos in Eq. J^l one write 



2K, 



■ cos 26* 



(9) 



is satisfied. This resonance converts a left-handed elec- 
tron neutrino into a right-handed (sterile) muon neutrino 
for the Dirac case and into a muon anti-neutrino in the 
Majorana case. (In principle there is another resonance 
possible for Dirac neutrinos, converting the chirality of 
the electron neutrino, but keeping its fiavor the same 
through a diagonal moment. But for the neutrinos that 
also go through the MSW resonance -as higher energy 
solar neutrinos do- this requires a very high neutron den- 
sity, Nn = 2Ne, which is not realized in the Sun). Clearly 
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TABLE I: The location of the MSW and spin-flavor pre- 
cession (SFP) resonances for Majorana neutrinos. Neutrino 
energies are given in MeV. The v/Rq value for the location 
of the resonances are shown. 



SFP 



MSW 



2.50 
3.35 
5.00 
8.00 
13.00 




0.05 
0.10 
0.15 
0.20 



0.07 
0.10 
0.13 
0.18 
0.22 



the resonances of (jH)) and Q that flip both the chirahty 
and the flavor of the electron neutrino produced in the 
nuclear reactions at the solar core take place at a higher 
electron density than the MSW resonance density. Even 
though their locations are different for the Dirac and Ma- 
jorana cases, neutrinos need to go through them first. To 
estimate its location we use the approximate expression 
for the solar electron density 



Neir) = 245exp(-10.54r/i?Q)iVA cm" 



(10) 



where Na is the Avogadro's number. For the solar neu- 
tron density we use a spline fit to the values given in 
Ref. IHl. Using the values of 5m^ = 8.2 x lO^^^ eV^ and 
tan^ 6 = 0.4, obtained from a global analysis of the solar 
neutrino and most recent KamLAND data ^2 we calcu- 
late the location of both the spin-fiavor and MSW reso- 
nances. Results for the Majorana neutrinos are given in 
Table I. One can see that since the magnetic field should 
be present at the location of the spin-flavor precession 
resonance only fields at and very near the core play a 
role. 

In this paper we present several approximate formulae 
for the electron-neutrino survival probability in several 
limiting cases. We first consider the case of a small mix- 
ing angle (non-adiabatic limit). In this case the SFP and 
MSW resonances are well separated. The derivation of 
the reduction of the electron neutrino survival probabil- 
ity is presented in detail in the appendix. Following Eq. 
(|54|l the electron neutrino oscillation probability in pres- 
ence of a magnetic field can be rewritten as 



X exp [—■7T{fiB)Sr] , 



P{Vel VeL.fJ'B = 0) 

(11) 



where the width of the SFP resonance (see Eq. (|50|l of the 
appendix) is given as 



5r 



K-K 



5m? cos 26 



(12) 



Clearly even for the rather large values of the magnetic 
field B ~ 10^ G and the magnetic moment of 10~^^ 
for a 10 MeV neutrino the width of the spin-fiavor res- 
onance would be very small, i.e. (^"^^ ~ 0.002. It is 
worth to reiterate that the value of the magnetic field at 



the close vicinity of the solar core is not restricted by he- 
lioseismology as our approximations are no longer valid 
for very large values of the magnetic field. 

It is not possible to find an expression similar to Eq. 
(|11|1 for larger mixing angles or the adiabatic limit in 
which case the SFP and MSW resonances are no longer 
well-separated. However it may be beneficial to inves- 
tigate the unrealistic, but pedagogically instructive limit 
where the neutron density vanishes. In this limit the SFP 
and MSW resonances overlap for any values of the neu- 
trino parameters and magnetic field. From Eq. (|24|l one 
can calculate 



,2 , -dfl) 



flBV2GFNn 



i^^Bf 
-- 



(13) 



In writing the above equation we assumed that the mag- 
netic field is a constant. In the limit Nn goes to zero the 
above equation becomes an equation for vi^^ only and 
can easily be solved using the semiclassical methods of 
Ref. For large initial electron densities we obtain 

Pil^e ^^^)^\^\ COS (1 - 2Phop) , (14) 

where the hopping probability is given by 

Phov = cxp(-7rfi), (15) 

with 



17 = - 



dr- 



2E 



[(C'(r)-2C(r) cos2a„ + l) 



{^,B) 



21 1/2 



(16) 



where Tq and ro are the turning points (zeros) of the 
integrand. In this expression we introduced the scaled 
density 



2V2GFNe{r 
5iin? / E 



(17) 



Using the Taylor series expansion of Eq. H16(l for small 
values of ^lB we can write 



exp ■ 



dr 



Phopi^^B ^ Q) = PhopitiB = Q) 
'2E ^C^(r)-2C(r) cos26l„ + l 



(18) 



For an exponential density, Ne{r) 



noe 



the inte- 



gral above can be calculated to give a hopping reduction 
factor 



exp 



TT {^iBf2E 
a 5m? 



(19) 



For a 10^ G magnetic field, a magnetic moment of 
lQ~^^^Bi and E ^ lOMeU, this hopping reduction fac- 
tor is very small, ^ 10^'^. For the near-adiabatic limit, 
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FIG. 1: Electron neutrino survival probability as a function 
of neutrino energy, in presence of both the MSW and the RSF 
resonances. In particular, a gaussian magnetic field profile is 
taken, with a magnetic field of 10^ G and a magnetic moment 
of 10~^^ hb- The location of the resonances are given in Table 
I. The results shown are indistinguishable from those due to 
the MSW resonance only, at least at the level of 10^. 



where we defined 

Q{t) = / dt'[cl){t') + K{t')]. (22) 
Jo 

Using Eq. (pn|l , Eq. (|^ can be easily integrated to yield 
n = 2^?L [ dtBl{t), (23) 

which can clearly be very large depending on the chosen 
average magnetic field profile. 

In conclusion, neutrinos emitted by the sun undergo a 
spin-flavor precession resonance and then an MSW res- 
onance. In this paper we have discussed the conditions 
to be met to encounter such resonances and derived ana- 
lytical formulae for the reduction of the electron survival 
probability due to solar magnetic fields. Our results show 
that the coupling of the neutrino magnetic moment to 
the solar magnetic field have small effects on the neu- 
trino fluxes. Such results indicate that future solar neu- 
trino measurements could not easily reach the level of 
precision to pinpoint alternative solutions to the solar 
neutrino deficit than the oscillation one confirmed by the 
recent KamLAND data. 



not only the hopping reduction factor is very small, but 
also the hopping probability itself is significantly reduced. 
Consequently one expects the change in the electron neu- 
trino survival probability due to a near-allowed value of 
the neutrino magnetic moment to be very small. We do 
not expect using a realistic, non-zero value of the neutron 
density to change this conclusion. For Nn ^ 0, however, 
an analytic expression does not exist. Therefore, for the 
case of the observed large mixing angle, we present results 
of a numerical calculation, obtained by solving directly 
Eq.(24). Figure 1 shows the electron survival probability 
as a function of the electron neutrino energy in presence 
of a magnetic field of 10^ G and a magnetic moment of 
10~^^/xb. In particular a gaussian profile has been taken 
for the magnetic field. Note that such results are identi- 
cal (at a level at least of 10^) to those obtained with the 
MSW effect only. This result is not changed if a width 
twice or five times larger than Eq.(12) is taken. 

It was recently argued in Ref. [i^ that random mag- 
netic fields can increase the reduction factor by en- 
hancing solar antineutrino flux. It is possible to show 
that this conclusion also follows from our approximate 
expressions. We assume a random magnetic field with 
the correlation function 

< B{r)B{r') >= B^{t)LS{r - r') (20) 

where Bq is the average value of the magnetic field and L 
is the correlation length. We rewrite the reduction factor 
given in Eq. (|53|l in the form 

n^2fi^ [ dtB{t)e"^^''> [ dt'B{t')e~'^^'\ (21) 
Jo Jo 
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Appendix: Reduction of the electron neutrino 
survival probability in the presence of magnetic 
fields 

In order to obtain the reduction of the electron neu- 
trino survival probability due to solar magnetic fields we 
implement the logarithmic perturbation theory (for a de- 
scription see the Appendix of Ref. ^3)- Here we show 
the derivation of the probability reduction in the case of 
Majorana neutrinos. A similar derivation can be given for 
the Dirac case. Note that the results presented here can 
also be used to explore the effect of the magnetic fields on 
the neutrino fluxes produced during the explosion of core- 
collapse supernovae. Since, as described above, there is 
no resonant production of i^i^'"^ in the Sun, we will set its 
amplitude to zero in Eq. (^) to obtain 

= A -0 (24) 

where we defined 

A = ^sin20, (25) 
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6m 
Is 



cos 26*, 



and 



Introducing 



and 



> 



(26) 



(27) 



(28) 



(29) 



Eq. can be rewritten as two coupled nonlinear equa- 
tions: 



I ; 

dt 

and 



A - 20Z 



(L) 



(30) 



dt 



= AiS-(«;+,^)z(«)-Az(^)z(^)-/iB[z(«)]2. (31) 



From Eqs. (jSOI) and (jSU it follows that 



'dt ( ' 



|^(L)|2^, (fl)|2 



= ^iB - z(«)) + A - . (32) 

Using the unitarity of the neutrino amplitudes (i.e. 1 + 



U(i)|2 + U(K)|2 



|i^fV)we obtain an exact expression 



for the electron neutrino amplitude at any location inside 
or outside the Sun: 



|z.W(T)|^ 



exp s * y 



X exp j dt 



^iB{t) (z(-^)*(t) - z(^)(i) 



it) 



(t) 



(33) 



Eq. H33() . which is an exact result, represents the for- 
mulation of the N-flavor (or antiflavor) propagation prob- 
lem in the SU {N) / SU {N - I) x U{1) coset space instead 
of the usual SU (N) one. To illustrate its utility in sepa- 
rating different contributions we define the perturbation 
parameter 
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fiBo 4EnBo 



5m? sin 29 ' 



(34) 



which we will take to be small. In Eq. I|34() Bq is the 
maximal value of the magnetic field. Wc write the mag- 
netic field in dimensionless form 



(3 = B/Bo. 



(35) 



Defining a new variable r = A<, Eqs. H30|l and H31() take 
the form 

^ 1 _ 2^z^^~> - [2(^)]2 - .g/3z(^)z(«), (36) 
ar A 



and 
I— — = gP 



dz^'^'> ^ K + (t) 



(37) 



These equations need to be solved with the initial condi- 
tion z^^^ = = z^^-* at t = 0. We consider a perturbative 
solution of the form 



and 



(L) (L) , (L) , 2 (L) , 



(38) 



(39) 



Clearly Zq^\ when substituted into Eq. H33|l gives the 
MSW solution and the terms proportional to various 
powers of g are the corrections due to the existence of 
the spin- flavor precession. 

The quantity Zq^'' satisfies the equation 
whereas the evolution of the zj^' is given by 



(40) 



,■^^1 - ^(L) _yJL) (L) 

dT ~ A 1 ° 1 ■ 



(41) 



Eq. (|41l) implies that that the quantity zj^-* is a con- 
stant times an exponential. The only way to satisfy the 
initial condition is to set this multiplicative constant to 
zero. Hence the lowest order correction to z^^'' is Z2^\ 
satisfying the equation 



. dz. 



(L) 



. -2jzf)-2z(-)zf)-/3z(-)zr. (42) 

Similarly, as it is expected on physical grounds, zjf''^ van- 
ishes. The lowest order correction to z^^^ is given by zj^^-* , 
which satisfies the equation 

1—^^13 —z\'~z'^>z\'. (43) 

The solution of Eq. 14311 is given by 

gz[^\T) = -ie^' dt'l4,(t')+nit')+Az'y(t')]] 

dtMi?(i)e'"^-^*'*''^^*'^'*'"^*'^+^'°''^*'^", (44) 

and the solution of Eq. H42I) is given by 

5zf)(r)=ze[^'r'^*'[^(*')+^^o''(*')ll (45) 

X rdi/.i?(i)4'^^(04''^(0e'~^'^o'*''*^*'^+''^°''**'^". 
Jo 
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When the magnetic field is set to zero Eq. (|33|l gives 
the neutrino survival probability to be 



f-T 

exp < I / dt 



Af4^)*(i)-4^)(i) 



(46) 



It is easy to see that this result, along with Eq. I|4()|l . 
represents a resonance. Introducing 



*(T) = exp 



-I I dt 
lo 



(Az(^)(t)+0(t) 



(47) 



one observes that \i^i^\T)\'^ = |4'(r)p. It follows from 
Eq. (|40|l that ^ satisfies the differential equation 



02 (t) 



Ad^ 



(48) 



The rate of change of the probability is maximized when 
the right hand side of Eq. (|48|l is an extremum, which is 
achieved when = 0. The width of this resonance (the 
MSW resonance) is F = ^/A, which corresponds to a 
spatial width of ASr = 2/T or 



5r 



2A 



MSW 



{d(l)/dt)^=o ' 



(49) 



A similar argument, applied to Eq. (|31|) gives the width 
of the spin-flavor resonance to be 



d / (p + K 



(50) 



As we mentioned earlier the spin-flavor precession res- 
onance takes place before the MSW resonance. In most 



cases the quantity is very small at the SEP resonance 
zone and can be neglected in Eqs. 14111 and (I42|l . In this 

approximation — and 

gz\ '(1 ) = —le Jo 



dt^iB{t)e^ 



-i /Jdt'[0(t')+K(t')]] 



(51) 



Substituting these in Eq. H33() we obtain 

|i.f)(r)p 



exp ■ 




T 

X exp \ i \ dt 



A(4'^*(t)-.r(t)) ' 



(52) 



P{Uel VeL.tiB 7^ 0) = P{Vel ^ VeL,f-l'B = 0) 



X exp 



dtfiB{t)e 



[t J*dt'[,P{t')+K.{t')]] 



(53) 



If the SEP resonance width is rather small one can cal- 
culate the integral in Eq. (|53() rather accurately in the 
stationary phase approximation. The stationary point is 
where the derivative of the argument of the exponent is 
zero, i.e. 0-1- k — 0, the SEP resonance point. One finally 
gets 



P{l^el VeL,tJ.B 7^ 0) = P{iyel VeLyfJ-B = 0) 



X exp 



|(i(0-f K)/fit|(0+„)^o 



(54) 
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